New explicit exact solutions for the Lienard equation 

and its applications 



o 

(N 



Gui-Qiong Xu * 



Department of Information Management, College of Management, 
' Shanghai University, Shanghai 200444, PR China 



in 



Abstract 

In this letter, new exact explicit solutions are obtained for the Lienard equation, and the ap- 
plications of the results to the generalized Pochhammer-Chree equation, the Kundu equation and 
the generalized long-short wave resonance equations are presented. 
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m ■ 1. Introduction 
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Nonlinear partial differential equations (NLPDEs) describe various nonlinear phenomena in 
natural and applied sciences such as fluid dynamics, plasma physics, solid state physics, optical 
^ ■ fibers, acoustics, mechanics, biology and mathematical finance. It is of significant importance to 

construct exact solutions of NLPDEs from both theoretical and practical points of view. Up to now, 
many powerful methods for solving NLPDEs have been proposed, such as the inverse scattering 
method[T], Backlund and Darboux transform [2] -[3], Hirota's bilinear method[4j, truncated painleve 
expansion method[5]-|10j. homogeneous balance method pT]. variational iteration method [12], ho- 
motopy perturbation method [KJ, tanh- function method [H], Jacobian elliptic function expansion 
method [15]- [19], Fan sub-equation method [20]- |22j. auxiliary equation method [23]- |25j. F-expansion 
method [26]- 1281 and so on. 



The last five methods mentioned above belong to a class of method called subsidiary ordinary 
differential equation method(sub-ODE method for short). The sub-ODE method which were often 
used the Riccati equation, Jacobian elliptic equation, projective Riccati equation, etc. In this letter, 
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we choose the Lienard equation 



a"{C) + ^a(0 + ma^iC) + na^{^) = 0, Imn 7^ 0, 



as the subsidiary ordinary differential equation. By means of some proper transformations, a 
number of NLPDEs with strong nonUnear terms can be reduced to Eq.([T|), thus seeking expUcit 
exact solutions of these nonlinear equations can be attributed to solve ([1]). Therefore, to search 
for exact solutions of the Lienard equation ([T|) is a very important job and it has attracted much 
attention. For example, Behera and Khare[29] has shown that the exact solution of Eq.([TJ can be 
expressed in terms of the Weierstrass function. Dey et al.[3n] investigated Eq.([T]) and established the 
exact solution of a one-parameter family of generalized Lienard equation with pth order nonlinearity 
by mapping it to the field equation of the (/>^-field theory. By means of different methods, Kong[31j. 
Zhang [52]-[55] and Feng|34]-[55] have given some explicit exact solitary wave solutions of Eq.([T]). 
In Refs.|32]-|36). Zhang and Feng derived three kinds of solitary wave solutions of Eq.([T]) as follows: 

If ^ < 0, m > 0, n < or / < 0, m < 0, n > 0, Eq.([T]) possesses the solitary wave solution. 



2,1' 



Sm^ — 16n/ 



sech 



-1 + 



V2, 



m 



V3 



Wnl 



sech 



If Z < 0, m > and 3m? — Wnl = 0, Eq.([T]) admits exact solutions. 



«2(e) 



-— (1 + tanh(\/^^~ 
m 



-— (1 - tanh(\/^e~ 
m 



(3) 



The various methods used in [29]- [36] are very useful and the applications of the solutions of 
the Lienard equation to some important NLPDEs are quite perfect. However, it is natural to ask 
whether Eq.([T]) can support other new exact solutions. The present letter is motivated by the 
desire to improve the work made in |31j-|36j by introducing more solutions of Eq.([T|) including all 
the solutions given in [31] -[3^ but also other formal solutions. 

The rest of this letter is organized as follows. In Section 2, we find some new exact solutions 
for the Lienard equation ([T]). In Section 3, we use these special solutions to solve the generalized 
Pochhammer-Chree equation, the Kundu equation and the generalized long-short wave resonance 
equations. And we conclude the letter in the last section. 
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2. New exact solutions of the Lienard equation 



Generally speaking, it is difficult to give the general solution of Eq.([T]). In what follows, we will 
consider some special cases. Based on Refs.[31]-[36], we can have the following solutions of Eq.([T]), 



ai±(0 

a2±(e) 

a3±(e) = ± 
a4±(0 = ± 



m + e^w? — 16n//3 cosh(2-v/^^) 
-M 



- 16n//3 > 0, / < 0, 



ra + ey^l6n//3 — sinh(2A/^^) 
— ( 1 + e tanh(\/^0 

(^1 + e coth(v^O 
-4/ 



, - 16n//3 < 0, / < 0, 



2/ 



, - 16nV3 = 0, m > 0, / < 0, n < 0, 
, - 16n//3 = 0, m > 0, / < 0, n < 0, 



- 16nV3 > 0, / > 0, 



(4a) 

(4b) 

(4c) 
(4d) 

(4e) 



e^rri^ - 16n//3 cos(2^/I ^) 
where e = ±1. It is easily seen that a3±(^) reproduces two solutions given in Eq.(l3]). There is a 
tiny symbolic error in the solution ai(^)( the coefficient of sech'^ \f—l£, in the numerator of fraction 
([2]) should be — 4-y^ J ^-i^ni )■ easily proved that the correct solution ai(^) and the solution 
oi-i-(^) with e = 1 are actually the same and only different in the form. And the other solutions 
o^2±(0) '^4±(0 ^-nd a^±{£,) are firstly reported here. 

To our best knowledge, the periodic wave solutions expressed in terms of Jacobian elliptic 
function to Eq.([T]) have not been considered in existed literature. Now we assume JacobiSN(^, r) = 
sn(^), JacobiCN(^, r) = cn(^) and JacobiDN((^, r) = dn(^), and r is the modulus of Jacobian elliptic 
functions(0 < r < 1). With the aid of symbolic computation software such as MAPLE, after direct 
computations, we find three kinds of elliptic periodic wave solutions of Eq.([T]) when the parameter 
coefficients /, m, n satisfy certain conditions. 



a6±(6 



a7±(0 



a8±(0 



3m / 



3m- / 

-— 1 + ecu 



^/3 



m 



n 



8n 



3m / , 

1 + edn 

8n I 



4r 

("v/S m 
4r ^/n 

-v/3 m 



3m2(5r2 - 1) 
^ = 2 — ^ m > 0, n < 0, (5a) 



, Sm'^Ur'^ + 1) 

/ = — -, m < 0, n > 0, (5b) 



64:12 



Z = ^ m < 0, n > 0. 5c 

64n ' ^ ^ 



To our knowledge, the solutions a6±(C)) 07±(^) and a8±(C) are firstly presented here. 



It is well known that there are many other Jacobian elliptic functions which can be generated 
by sn(^), cn(^) and dn(^). For the sake of simplicity, the solutions in terms of ns(^), nd(^), nc(^), 
sc(^), cs(^), sd(^), ds(^), cd(^), dc(^) are not considered here. 



3. Applications 



Example 1. The generalized Pochhammer-Chree (PC) equation can be written as 



Utt - Uttx 



(ai u + a^v? + a^u^ 



0, 



(6) 



which describes the propagation of longitudinal deformation waves in an elastic rod [37]. Zhang [32] 
and Feng[3S] have given some explicit solitary wave solutions of Eq.Q by means of the method of 
solving algebraic equations. Li and Zhang[38] studied the bifurcation problem of travelling wave 
solutions for Eq.([6]) by using the bifurcation theory of planar dynamical systems. 

In order to solve Eq.®, its solutions may be supposed as: 



u{x,t) = u(^), ^ = X — vt, 
where v is a real constant. Substituting ansatz d?]) into Eq.® yields, 

v^u"{^) - v^u^'^\0 - {aiu + agii^ + a^u^)^^ = 0, 
Integrating Eq.® twice and setting the integration constant to zero, we obtain 



(7) 



(8) 



(9) 



Up to now, by means of the ansatz ([7]), we reduce the generalized PC equation ([6]) to the Lienard 

2 

equation ([1]) for the case I = , m = ^ and n = p-. Substituting the solutions ([Ia|) - ([Ie|) and 
the solutions (f5a|l - (f5c]) of Eq.([T]) into ([7]), we can obtain a series of exact travelling wave solutions 
to Eq.® (where ei = ±1 and €2 = ±1). 

When t;^ — ai > and 3a| — 1605(01 — v'^) > 0, Eq.Q has bell-shape solitary wave solution. 



ui±{x,t) 



4{v^ - 01) 



03 + ei Y^ol — 1605(01 — v'^)/3 cosh( 



2^/l 



ai 







When — oi > and 3a| — 1605(01 — f^) < 0, Eq.Q has the singular solitary wave solution, 



U2±{x,t) 



01 



03 + eiY^16a5(oi — f2)/3 — O3 sinh( 







A 



When as > 0, 05 < 0, — ai > and 803 — 1605(01 — v'^) = 0, Eq.(l6]) has two kink-shape 
sohtary wave solutions, 



U4±{x,t) = ± 



2{v^ - ai) 
03 



1 + ei tanh 



— ai 



When — ai < and 803 — 1605(01 — u^) > 0, Eq.([U]) has the trigonometric function solution, 

4(t;2 -01) 



U5±ix,t) = ± 



03 + eiy^Og — 1605(01 — v'^)/3 cos 



.2^/al-v^ 







When 05 < and 03 > 0, Eq.Q has the Jacobian sine function solution, 



/ , 1 



303 / / \/303 

1 + ei sn I -; ^= 4 



2 05 \ ' ^ \4rv \/^a5 

where = £2 \/ 0.5 (64 osr^oi — 15 os^r^ + 3 032)7(8 r 05). 

When 05 > and 03 < 0, Eq.([6]) has two periodic wave solutions. One is 

\/3a3 



U7±{x,t) = ±- 



303 



1 + ei cn 



205 y \ 4ri;^/a5 

where v = €2 \/ 0-5 (64a5r2ai — lla^'^r'^ — 3a^'^)/{8ra5). And another one is 



us±ix,t) = ±- 



' l + eidn(^e 



205 



where v = e2 \/ 05 (64 0501 — 3 os^r^ — 12 03^)/ (8 05). 

Among the above solutions, only ■ui±(x,t) with ei = 1 and U3±{x,t) reproduce the results given 
in Refs.[32]-[3S], and the other solutions have not been found before. 
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Example 2. Next we consider the Kundu equation, 

iut + Uxx + fi\u\^u + (5|ti|^ii + ia{\u\^u)x + i s {\u\^)xU = 0, 



where /3,(5, a, s are real constants. Eq. fjlOp was derived by Kundu[39] in the study of integrability 
and it is an important special case of the generalized complex Ginzburg-Laudau equation [40]. 
Meanwhile, Eq. (|10p and its special cases arise in various physical and mechanical applications, 
such as plasma physics, nonlinear fluid mechanics, nonlinear optics and quantum physics. Fengji 



derived the explicit exact solitary wave solutions of Eq. ijlOp by using the algebraic curve method. 
Zhang et a/.^41j studied the orbital stability of solitary waves for Eq. (fTO|) by means of spectral 
analysis. 



Assume that Eq. ijlOp has solutions of the the form 

u{x,t) = 0(e)e^('^(«)-'^*\ ^ = X - vt, 



(11) 



where u and v are constants to be determined. Substituting Eq. (fTT]) into Eq. (fTO|) and then sepa- 
rating the real part and imaginary part yields, 

(CO + vi;'{o)m + <t>"{o - m^y^o - «0'(ov''(o + i3<pHo + = o, (12a) 



(12b) 



Letting 

i^'iC) = A + B<t>\0- (13) 

Substituting Eq.l^ into Eq.jUbl) and setting the coefficients of (/>'(C), (j)'^{04>'{0 to zero, we have 
A = v/2, B = -{3a + 2s) /A. Then Eq.^ becomes, 

i'r^\ ^ 3a + 2s 2/^\ /I A\ 

^ (6 = 2 1 (0- (14) 

Substituting Eq. (jl4p into Eq. (jl2ap yields the Lienard equation of the form, 

f{.€) +im +m<l,\i) +n<l,\i) = 0, 



(15) 



where l,m,n are given by 



^ av , (a - 2s) (3a + 2s) 
I = UJ + -, m = f3-—, n = 5 + ^ '-^ 



By the transformations pT|) and the exact solutions of Eq. (|TO|) can be obtained by using 
the solutions of Eq.dl]) given in Section 2. In the following solutions, V'(C) is given by Eq. ()14p . 

Ai = (2/3 - au)2 - (4a; + v'^){165 + {a - 2s)(3a + 2s))/3. 

When Ai > 0, and f + 4w < 0, Eq. ()10p has the solitary wave solution, 
ui{x,t) = (t>{x- ^;t)e^W^-^*)-'^*), 



-2(4w + t>2) 



2/3- av + e^/A^ cosh(Y^-(u2 + 4^;) ^) 
When Ai < 0, and f ^ + 4w < 0, Eq. (fTO|) has the singular solitary wave solution. 



2(3 -av + eV-Ai sinh(Y^-(?;2 + 4a;) ^) 



When v'^ + Auj < and av — 2/3 < 0, Eq.([TO|) has two kink-shape sohtary wave solutions, 
U3{x,t) = (f)ix- ^;t)e^(V'(^-^*)-'^*), 



m = ± 



a v - 2/3 

Ui{x,t) = (j){x- t,t)e^(V'(^-''*)-'^*), 



l + .tanh(V3?^^) 



m = ± 



aw - 2/3 



i + ecoth(vE(l^e) 



where uj is determined by Ai = 0. 

When Ai > and v'^ + Aui > 0, Eg. ([TO]) has the periodic solution of trigonometric function, 

-2{4uj + 



2/3 - + e^AT cos(\/f ^ + 4a; ^) 



When 4sa + 4s^ — 3a^ — 16(5 > 0, 2f3 — av > 0, Eq. pO|) has the Jacobian elliptic sine function 
solution. 



<^(6 = ± 



3(2/3 - aw) 



4sa + 4s2 - 3a2 - 16(5 



1 + esn( 



V^{2/3-av) 



2r V4:sa + 4s2 - Sa^ - 165 







where w is determined by r'^{v'^ + 4w)(16(5 + (3a + 2s)(a - 2s)) - 3(/3 - va/2)'^{5r'^ - 1) = 0. 

When 4sa + 4s^ — 3a^ — 165 < 0, 2/3 — av < 0, Eq.pU]) has two Jacobian elliptic function 
solutions. One is 



U7{x,t) = (pix- t;t)e*(^(^-^*)-'^*), 
3(2/3 -aw) 



m = ± 



4sa + 4s2 - 3a2 - 166 



1 + ecn( 



\/3(2/3 



aw 



2r V3a2 + 16(5 - 4sa - 4s2 



6 



where u is determined by {v^ + 4a;)(16(5 + (3 a + 2 s)(a - 2 s)) - 3 (/3 - wa/2)2(4 + 1) = 0. And 
another one is 

U8{x,t) = (pix- t,t)e*(V'(^-''*)-'^*), 
3(2/3 - aw) 



</'(0 = ± 



4sa + 4s2 - 3a2 - 16(5 



1 + edn( 



^/3(2/3 - aw) 



2 V3a2 + 16(5 - 4sa - 4s2 







where oj is determined by {v^ + 4u}){166 + (3a + 2s)(a - 2s)) - 3 (/3 - wa/2)^(H + 4) = 0. 
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The solutions ui{x, t) with e = 1, ^3(2;, t), U4{x, t) are same as the results reported in [34]. Other 
solutions have not been reported in [34J. In addition, the Kundu equation (jlOp contains several 
important nonlinear models when taking different choices for the parameters a, /?, 6 and s. For 
example, if s = 0, Eq. pOj) reduces to the derivative Schrodinger equation [39] 

iut + Uxx + /3|wpn + 5|n|^ti + ia{\u\'^u)x = 0; (16) 

if 5 = 20-2, a = -2a, s = Aa, then Eq.(^ becomes the Gerdjikov-Ivanov equation [32], 

iut + Uxx + f3\u\'^u + 2a'^ \u\^u + liau^Ux = 0. (17) 

Obviously, the explicit exact solutions of Eq. (jl6p and Eq. (jl7p can be derived from the above 
solutions. 

Example 3. Finally we consider the generalized long-short wave resonance equations with strong 
nonlinear term, 

iSt + Sxx = aLS + j\S\^S + S\S\^S, 

(18) 

Lt + P\S\l = 0, 

where S is the envelope of the short wave, and L is the amplitude of the long wave and is real. The 
parameters a,f3,j and 5 are arbitrary real constants. Recently, Shang[l3] obtained several kinds 
of explicit exact solutions of Eq. (fT8|) . 

In order to seek the exact solutions of Eq. (|18p . we introduce the following transformation, 

S{x,t) = (^(x,t)e*('=^+'^*+«o), (19) 

where <j){x,t) is a real- valued function, and k and w are constants to be determined, ^0 is an 
arbitrary constant. Substituting Eq. (|19p into Eq. (|18p and then separating the real and imaginary 
parts yields, 

(pxx - i^ + k'^)(l) - aL(j) - 'jcp^ - 6(1)^ = 0, (20a) 

ct>t + 2k(l)x = 0, (20b) 

Lt + 2p(t>(t>x = 0. (20c) 

In view of Eq. ()20bp we suppose 

cP{x,t) = cP{0=cP{x-2kt + Ci), (21) 
where ^1 is an arbitrary constant. Therefore we also assume 

L{x, t) = V(0 = Hx-2kt + 6). (22) 



Substituting Eq.dH]) into Eq.gUc^ yields, 



where C is an integration constant. 

Substituting Eqs.(l2l])-([23]) into Eq.i^^, we have, 

+ im + m^'iO + n<t>\i) = 0, 

where the parameters Z, m, n are given by 

I = - {uj + k'^ + aC), 'fn = -{1 + 1^)1 n = 



(23) 



(24) 



(25) 



Similar to Example 1, by means of the transformations (I19p . (|2ip -(f23]). we can also reduce 
the generalized long-short wave resonance equations (llSp to the Lienard equation ([T]). Together 
with Eq. (j2ip and Eq. ()23p . substituting the solutions of the Lienard equation given in Section 2 into 



Eq. (|19p and Eq. (j22p yields abundant periodic wave solutions of the generalized long-short wave 

a/3 

resonance equations (fTSll . In the following eight sets of solutions, A2 = (7 H -)'^ — 16(5 (a; + /c^ + 

2 k 

aC)/3, e = ±1, and ^ = x — 2A;t + ^i with k being nonzero arbitrary constant. 

When A2 > and w + A;^ + q C > 0, Eqs. (fT8]l has a set of bell-shape solitary wave solutions, 

4/3(w + A:2 + aC) 



Li(x,t) 



Si{x, t 



2/c 7 - a /3 + 2A: e cosh(2Vw + A;^ + a CO 
4 (u; + fc^ + g C) 



-(7 + ttt) + e^A^ cosh(2Vu; + P + QCe) 



J, {kx+ujt+io) 



When A2 < 0, w + /c^ + Q C > Eqs. ffTSl) has a set of singular solitary wave solutions, 



L2{x,t) 
S2ix,t) 



4:/3{uj + P + aC) 



2k-f - a/3 + 2ke V^A^ sinh(2Vw + A;2 + a 

4{LU + P + aC) 

a /3, 



+ c, 



-(7 + ^) + sinh(2Vu; + A:2 + aC'0 



When A2 = 0, a; + + aC > 0, and 2^7 + a/3 < 0, Eqs. ffTS]) has two sets of kink-shape solitary 
wave solutions, 

= _ 2/3(a; + fc2 + a^) (1 + e tanhWu: + k'^ + aC ()) + C, 
2A;7 + ap \ J 



S3{x,t) 



ikjcj + e + aC) 
2^7 + a/3 



1 + e tanh(\/w + P + aC^ 



Q 



L4{x,t) = - 

S4{x, t) = ± 



2k^ + ap \ / 



4:k{uj + k'^ + aC) 



1 + e coth(\/w + A;2 + aCO 



J (fex+cjt+fo) 



2^7 + a/3 

When A2 > and uj + k'^ + aC < 0, Eqs. (fT8|) has a set of trigonometric function solutions, 

A(3{cj + k'^ + aC) 



2k-i-ap + 2k e ^/A^ cos(2v^-(w + k'^ + aC)C) 

4 (u; + fc^ + g C) 

a/3, 



-(7 + ^) + ^Va; cos(2V-(c^ + A;2 + aC)0 



When (5 > and A;(a/3 + 2A;7) < 0, Eqs. psp has a set of Jacobian ehiptic sine function solutions, 
3/3(a/3 + 2A;7) 



LQ{x,t) ■- 
SQ{x,t) = ± 



32A;2 6 \ 

3(a/3 + 2fc7) 
16k 5 



x/3(a/3 + 2fc7) ^ „ 
1 + esn( — ^) + C, 



1 + esn( 



\/3 (a/3 + 2/27) 



a/3. 



where w is determined by 6Ar'^6 {lo + k'^ + a C) — 3{5r'^ — 1) (7 H — )2 = 0. 

2k 

When J < 0, k{a/3 + 2A;7) > 0, Eqs.([T8]) has two sets of Jacobian elliptic function solutions. 
One is 



L7{x,t) = - 
S7{x,t) =± 



3/3 (a/3 + 2A: 7) / 
32F5 1 

3(a/3 + 2/c7) 
Wk5 



^/3(a/3 + 2fc7) ^ „ 
1 + ecn( ; — = 4) + C, 



8k' 



J (kx+u) i+5o) 



a/3 

where w is determined by GAr'^S (w + fc^ + a C) — 3(4r2 + 1) (7 H -)2 = 0. And another one is 



Lsix,t) = - 
S8ix,t) = ± 



3/3 (a/3 + 2/c 7) / 
32k^6 1 

3{a/3 + 2kj) 
16k 6 



2k 



1 + ,d„(-^iM±|Mj)l + c. 



1 + edn(- 



V3{al3 + 2kj) 
8k7^ 



3? (fcx+cjt+^o) 



a/3 

where w is determined by 645 {to + k'^ + aC) — 3(4 + r^) (7 -| — )2 = 0. 

With the aid of Maple, we have checked all solutions by putting them back into the original 
Equation. 



in 



4. Conclusions 

The Lienard equation is used to describe fluid-mechanical and nonhnear elastic mechanical 
phenomena. Moreover, a number of NLPDEs with strong nonlinear terms can be reduced to the 
Lienard equation by some proper transformations. Therefore, To search for new special solutions of 
the Lienard equation is a very important job. In this letter, we obtain eight kinds of explicit exact 
solutions for the Lienard equation, which include solitary wave solutions, periodic wave solutions 
in terms of trigonometric function and Jacobian elliptic function. By means of these solutions, we 
obtain a variety of explicit exact solutions for the generalized PC equation, the Kundu equation and 
the generalized long-short wave resonance equations. These solutions may be important explain 
some physical phenomena. The method presented here is also applicable to solve other nonlinear 
equations with strong nonlinear terms. For example, the Ablowitz eg nation |44j. 

iutt = Uxx — '^iu^Ux + 8|n|^tt; 

the third-order generalized NLS equation(also called RKL model) [45] . 

iuz + utt + 2\u\^u + iauttt + 'i'li{\u\^u)t + il{\uf'u)t + 5\u\'^u = 0, 

and the nonlinear equations which were considered in Ref . |32] and Ref. |34j . 
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